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Abstract 

The cross section for the thermal neutron capture by the deuteron is calculated with pionless 
<^> Effective Field Theory (EFT). No new Three-Nucleon forces are needed up to next-to-next-to- 
leading order in order to achieve cut-off independent results, besides those fixed by the triton 
^ binding energy and Nd scattering length in the triton channel. The cross-section is accurately 
Q determined to be atot = [0.503 ± 0.003]m6. At zero energies, the magnetic Ml-transition gives 
^ the dominant contribution and is calculated up to next-to-next-to- leading order (N 2 LO). Close 
agreement between the available experimental data and the calculated cross section is reached. 
We demonstrate convergence and cutoff independence order by order in the low-energy expansion. 
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I. INTRODUCTION 



The study of the three-body nuclear system involving neutron radiative capture by 
deuteron has been investigated in theoretical and experimental works over the past years. 
The experimental result of this process has most accurately been measured by Jurney, 
et.al. Q. The value of 0.508±0.015(mb) for the cross section was resulted for 2200 m/sec 
neutrons. 

Rapid progress has been made in the theoretical study of the Nd — > 3 i?7 reaction such 
as the p-d and n-d radiative capture. At such energies a magnetic dipole(Mi)transition is 
almost entirely participated. These reactions were studied in plane wave (Born) approxi- 
mation by Friar et al. [2j . In these investigations the authors employed their configuration- 
space Faddeev calculations of the helium wave function, with inclusion of three-body forces 
and pion exchange currents. More recently a rather detailed investigation of such pro- 
cesses has been performed by Viviani et al. [!, 19]. In their calculations the quite accurate 
three-nucleon bound- and continuum states were obtained in the variational pair-correlated 
hyperspherical method from a realistic Hamiltonian model with two- nucleon and three- 
nucleon interactions. 

They obtained in Ref. 0] the cross section from Argonne v 14 two-nucleon and Urbana 
VIII three-nucleon interactions(AV14/UVIII), also from Argonne Vi$ two-nucleon and Ur- 
bana IX three-nucleon interactions(AV18/UIX) and including A admixtures. Cross section 
values were found 0.600 (mb) and 0.578 (mb) which overestimate the experimental value by 
18% and 14% value, respectively, see table 2. It should be noted, however, that the explicit, 
non-perturbative inclusion of A-isobar degrees of freedom in the nuclear wave function are 
found to be in significantly better agreement with experiment than those obtained from per- 
turbative (Apt) estimates. This shows that these results for this very- low energy observable 
are sensitive to details of the short-range part of the interaction, recent calculation using 
manifestly gauge-invariant currents reduced the spread but the result including three- 
body currents, 0.558 mb, still over-predicts the cross-section by 10%. Model-dependent 
currents associated with the A(1232) were identified as source of the discrepancy. Thus, 
the question remains how such details of short-range Physics can so severely influence a 
very-long-range reaction with maximal energies of less than 10 MeV. 

During the last few years, nuclear Effective Field Theor y(E FT) has been applied to two-, 
three-, and four- nucleon systems, see e.g. 0, H, @, 0, @, BUS]- The pionless Effective Field 
Theory would be an ideal tool to calculate low-energy cross sections in a model-independent 
way and to possibly reduce the theoretical errors by a systematic, model-independent cal- 
culation with an a-priori estimate of the theoretical uncertainties. An example of a precise 
calculation is the reaction np — > jd, which is relevant to big-bang nucleosynthesis(BBN). 
The cross section for this process was computed to 1% error for center of mass energies 



E < lMev 11, 12, 13 



We have suggested a method for computation of neutron-deuteron radiative capture 
for extremely low energy( 20 < E < 200 Kev )with pionless EFT 15j, where with this 
formalism, we can estimate errors in a perturbative expansion up to N 2 LO within a few 



percent of the ENDF values [16 



The purpose of the present paper is to study the cross section for radiative capture of 



2 



neutrons by deuterons nd — > 7 3 if at zero energies with pionless EFT. At these energies, the 
magnetic Mi-transition gives the dominant contribution. The Mi amplitude is calculated 
up to next-to-next-to-leading order (N 2 LO) with insertion of three body force. Results show 
less than 1% deviation from the available experimental data at zero energy (0.0253 eV). 

This article is organized as follows. In the next section, a brief description of the formal- 
ism and its input for total cross section of the neutron-deuteron radiative capture will be 
presented. We discuss the theoretical errors, tabulation of the calculated cross section in 
comparison with the other theoretical approaches and the most recent data [l[ in section [TTTl 
Finally, Summary and conclusions follow in Section IIVI 



II. NEUTRON-DEUTERON SCATTERING IN TRITON CHANNEL AND RA- 
DIATIVE CAPTURE 



The 2 Si channel to which 3 He and 3 H belong is qualitatively different from the other 
three-nucleon channels because all three nucleons can occupy the same points in space. 
Consequently, 2 Si/ 2 describes the preferred mode for nd — > 3 Hj and pd — > 3 Hej. The 
three-nucleon Lagrangean is well-known and will not be repeated here, see e.g. [14,18] for 
details. 

The derivation of the integral equation describing neutron-deuteron scattering has also 



been discussed before, see e.g. [7|,ll8|]. We present here only the results. The integral equation 



is solved numerically by imposing a cut-off A. In that case, a unique solution exists in the 
2 5i/2-cliannel for each A and vanishing three-body force, but no unique limit as A — > oo. 
As long-distance phenomena must however be insensitive to details of the short- distance 
physics (and in particular of the regulator chosen), Bedaque et al. jl, 0, 0, IS] showed that 
the system must be stabilized by a three-body force 



UiE; A) = £ ± "-(A) (^) " = ^ + ^ [ME + tf) + 

n=0 ^ ' 



(1) 



which absorbs all dependence on the cut-off as A — > oo. It is analytical in E and can be 
obtained from a three-body Lagrangean, employing a three-nucleon auxiliary field analogous 
to the treatment of the two-nucleon channels [14| . Contrary to the terms without derivatives, 
there are different, inequivalent three-body force terms with two derivatives, but only one 
of them, H 2 , is enhanced over its naive dimensional estimate, mandating its inclusion at 
N 2 LO [14,20]. Neutron-deuteron scattering amplitude including the new term generated by 
the two- derivative three-body force is shown schematically in Fig.l. Two amplitudes get 
mixed: t s describes the d t + N — > d s + N process, and t t describes the d t + N d t + N 
process, where d s (d t ) is an auxiliary field of two nucleons in a relative singlet-S (triplet-S) 
wave. 
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FIG. 1: The Faddeev equation for iVd-scattering. Thick solid line is propagator of the two inter- 
mediate auxiliary fields D s and D t , denoted by T>; K,: propagator of the exchanged nucleon; 7i: 
three-body force. 



A 



t s (p,k) =\[3K(p,k) + 2H(E,A)) + ± J dqq 2 [V a (q)[K(p,q) + 2H(E, A)]t s (q) 

o 

+V t (q) [3Kfaq) + 2H(E,A))tt(q)] 



A 



t t (p, k) = \ [Kip, k) + 2H(E, A)] + i- J dq q 2 [V t (q) [IC(p, q) + 2H(E, A)] t t (q) 

o 

+V s (q) [3/C(p, q) + 2H(E, A)] t s (q)) (2) 

2 

where V Stt (q) = T> Sit (E — j^,q) are the propagators of the auxiliary fields d Sjt , and /C the 
propagator of the exchanged nucleon, projected into the S-wave. For the spin-triplet S-wave 
channel, one determines the two-nucleon interaction up to N 2 LO by the deuteron binding 
momentum j t = 45.7025 MeV and effective range p t = 1.764 fm. Because there is no real 
bound state in the spin singlet channel of the two-nucleon system, its free parameters are 
better determined by the scattering length a s = 1/js — —23.714 fm and the effective range 
r s = 2.73 fm at zero momentum. 

The neutron-deuteron J = 1/2 phase shifts 6 is determined by the on-shell amplitude 
t t (k, k), multiplied with the wave function renormalisation 

T {k) = Z<t (*, k) ^—±—. (3, 

At thermal energies, the reaction proceeds through 5*-wave capture predominantly via 

a magnetic dipole transition, M^ SJ , where L=0, 5=1/2,3/2 and i=l. To obtain the spin 

structure, which corresponds to a definite value of J for the entrance channel, it is necessary 

- 1+ 3 + 

to build special linear combinations of products DN and axDN, with J p = - or J p = — , 

and D the deuteron spin-one field, see [15] for details. 



h/2 



{iD + ax D)N and [2iD -ax D)N 



1+ 3- 
For both possible magnetic dipole transitions with J p = — (amplitude g\) and J p = — 
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(amplitude g 3 ) we can write: 



#i : t ] {iD ■ e* x k + a x D ■ e* x fc)iV, 

# 3 : £ f (iD • e* x k + a x £> • e* x fc)7V . (4) 

The contribution of the electric transition E^ SJ for energies of less than 60 KeV to the 
total cross section is very small. Therefore, the electric quadrupole transition e®( 3 / 2 ^ 3 / 2 ^ 
from the initial quartet state will not be considered at thermal energies. The M\ amplitude 
receives contributions from the magnetic moments of the nucleon and dibaryon operators 
coupling to the magnetic field, which are described by the Lagrange density 

C B = ^-N\k, + hr 3 )a.B + e ^ d/d s3 B j + H.C . (5) 

ZM N M N Vr( S °)H Sl > 

where k = l/2(k p + k n ) = 0.4399 and k\ = l/2(k p — k n ) = 2.35294 are the isoscalar and 
isovector nucleon magnetic moment in nuclear magnetons, respectively. The NLO-coefficient 



Li is fixed at its leading non-vanishing order to the thermal cross section [11 




a-B 

FIG. 2: Some diagrams for adding photon-interaction to the Faddeev equation up to N 2 LO. Wavy 
line shows photon and small circles show magnetic photon interaction. For L\ vertices, see eq.(7); 
f^2 : three- body force, see eq.(3). Remaining notation as in Fig. [TJ 



The radiative capture cross section nd — > 3 H'j at very low energy is given by 

.3 



a 



9 v rel AM N 



E0xf 5 T], (6) 



r2 

' iLSJ 



where 

X\ SJ = —V^ X f SJ , (7) 

PUN 

with x stands for either E or M and hn is in nuclear magneton and p is momentum of the 
incident neutron in the center of mass. 

We now turn to the Faddeev integral equation to be used in the Mi calculation. We solve 
the Faddeev equation for nd-scattering and also for the triton bound state to some order (e.g. 
LO), then we take these Faddeev amplitudes and sandwich the photon-interactions with 
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nucleons between them when the photon kernel is expanded to the same order. This process 
will be done separately for NLO and N 2 LO. Finally the wave function renormalization in 
each order will be done. 

The diagrams in Fig. [2] represent contributions of electromagnetic interaction with nu- 
cleon, deuteron, four-nucleon-magnetic-photon operator described by a coupling between 
the 3 5*i-dibaryon and 1 5'o-dibaryon and a magnetic photon. As mentioned in the intro- 
duction, in another paper [15], we have presented detailed schematic of these diagrams in 
neutron-deuteron radiative capture for (20 < E < 200 keV) up to N 2 LO. 

The last diagrams in Fig. [2] with insertion of a photon to the N 2 LO three-nucleon force 
H 2 vertex is not Mi and we know that Mi contribution is the dominant contribution at 
very low energy and especially for zero energy. Its contribution should therefore be very 
tiny. Because the leading three-nucleon force H has no derivatives, it is not affected by the 
minimal substitution p — > p — e A But the parameter H 2 is the strength of the three-nucleon 
interaction with two derivatives. Naturally for the energy range near zero momentum, 
insertion of photon to H 2 vertices for momentum p ~ 0.025eV and Mi transition, could be 
neglected. H 2 is necessary in neutron-deuteron scattering to improve cut-off independence 
but is defined such that it does not contribute at zero momentum. Contributions of a photon 
coupling to H 2 are however indeed negligible at zero energy. 

III. NEUTRON-DEUTERON RADIATIVE CAPTURE RESULTS AT ZERO EN- 
ERGY 

We numerically solved the Faddeev integral equation up to N 2 LO. We used he = 
197.327 MeVfm, a nucleon mass of M = 938.918 MeV, for the NN triplet channel a 
deuteron binding energy (momentum) of B = 2.225 MeV (7^ = 45.7066 MeV), a residue of 
Zd = 1.690(3), for the NN singlet channel an 1 So scattering length of a s = —23.714 fm , 
Li ~ —4.5 fm by fixing at its leading non-vanishing order by the thermal cross section. 

As in Ref. [20], we can determine which three-body forces are required at any given order, 
and how they depend on the cutoff. 

Low-energy observables must be insensitive to the cut-off, namely to any details of short- 
distance physics in the region above the break-down scale of the pion-less EFT, set approx- 
imately by the pion-mass. It was found in Ref. [20] that no additional three-nucleon forces 
are necessary to render a renormalisable amplitude at N 2 LO in this process, besides those 
needed already in nucleon-deuteron scattering: H and H 2 . At N 2 LO , where we saw that 
H 2 is required, we checked this by varying the cut-off between 150 and 500 MeV. This is a 
reasonable estimate of the errors of our calculation due to higher-order effects. As seen in 
Fig. 3, in the thermal energy range the cutoff variation is very small and decreases steadily 
as we increase the order of the calculation and it is of the order of (k/A) n , (7/A) n , where 
n is the order of the calculation and A = 150 MeV is the smallest cutoff used ( see Table |T] 
and Fig. [3]). Also, errors due to cutoff variation is decreasing when the order of calculation 
is increased up to N 2 LO. 

We determined the two-nucleon parameters from the deuteron binding energy, triplet 
effective range (defined by an expansion around the deuteron pole, not at zero momentum), 
the singlet scattering length, effective range (defined by expanding at zero momentum), and 
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FIG. 3: Curve of the cutoff variation of cross section up to N 2 LO is shown between A = 150 MeV 
and A = 500 MeV. The short dashed, long dashed and solid line correspond to LO, NLO and 
N 2 LO , respectively. 

TABLE I: Results for the cutoff variation of cross section up to N 2 LO is shown between A = 150 
MeV and A = 500 MeV. 



E(10~ 8 MeV) 


LO 


NLO 


N 2 LO 


1 


0.0006 


0.00005 


0.0000002 


2 


0.0010 


0.00040 


0.0000005 


2.65 


0.0012 


0.00060 


0.0000090 


3 


0.0014 


0.00084 


0.0000150 


10 


0.0020 


0.00131 


0.0000600 



two body capture process(obtained with comparison between experimental data and theo- 
retical results for np — > dj process at zero energy [12]). We fix the three-body parameters as 
follows: because we defined H 2 such that it does not contribute at zero momentum scatter- 
ing, one can first determine H from the 2 Si scattering length = (0.65 ±0.04) fm [17]. At 

LO and NLO, this is the only three-body force. At N 2 LO, H2 is required. It is determined 
by the triton binding energy B% = 8.48 MeV. Finally, we solve by insertion of the potential 
at a given order in the integral equation and iteration of kernel. 

The cross section for neutron-deuteron radiative capture as function of the center-of-mass 
energy up to N 2 LO is shown in Fig. |H We also show single point that shows the available 
experimental results for this cross section at 0.025 eV [lj]. 

Table [Til shows Comparison between results of different models-dependent, model- 
independent EFT and experiment, for neutron radiative capture by deuteron up to N 2 LO, 
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FIG. 4: The cross section for neutron radiative capture by deuteron as function of the center-of- 
mass kinetic energy E in MeV. The short dashed, long dashed and solid line correspond to the 
contribution of M\ capture cross section up to LO, NLO and N 2 LO , respectively. Single point 
shows experimental results for this cross section at 0.025 eV 

TABLE II: Comparison between different theoretical results for Neutron radiative capture by 
deuteron at zero energy (0.0253 ev). Last row shows our EFT result. The last line quotes deviation 
between data [1] and theory, if it is larger than the theoretical or experimental uncertainty. 



Theory 



cr(mb) 



deviation from exp. 



AV14/VIII (IA+MI+MD) [3] 
AV18/IX (IA+MI+MD) [3] 
AV14/VIII(IA+MI+MD+A PT ) [3] 
AV18/IX(IA+MI+MD+A PT ) [3] 
AV14/VIII(IA+MI+MD+A) [3] 
AV18/IX(IA+MI+MD+A) [3j 
AV18/IX (gauge inv.) [19] 
AV18/IX (gauge inv. + 3N-current) [19] 



0.509 
0.489 
0.658 
0.631 
0.600 
0.578 
0.523 
0.556 
0.485 
0.496 



4% 
29% 
24% 
18% 
14% 
3% 



EFT(LO) 
EFT(NLO) 
EFT(N 2 LO) 
Experiment [1] 



0.503 ± 0.003 
0.508 + 0.015 



5% 




The calculations by Viviani et al. 0, 19[ shows sensitivity to 
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short-range physics namely to details of including the physics of the Delta and pion-exchange 
currents. The calculation of Ref. [19] with manifestly gauge- invariant current operators is 
quite sensitive to including meson-exchange three-nucleon currents. One might therefore 
have been tempted to conclude that a new three-nucleon force is also needed in the pion- 
less EFT. As shown above, this is not the case: There are no new three-nucleon forces 
besides those already fixed in nd scattering at the same order. The contribution from the 
photon coupling to a three-nucleon force is negligible in our calculation. As our result is 
model-independent and universal, any model with the same input must - within the ac- 
curacy of our calculation - lead to the same result. Our inputs are the first two terms of 
the effective-range expansion in the singlet- and triplet-S wave of NN scattering, the proton 
and neutron magnetic moments, the triton binding energy and nd scattering length in the 
doublet-S-wave, and finally the thermal cross section of the reaction np — > (determining 
L\). More work is needed to understand why the potential- model calculations [3,19] have 
the same input but do not seem to reproduce the same result. 

Addressing convergence of the EFT calculation, we notice that the contributions which 
are characterized as higher-order in the power-counting are indeed small: The LO result is 
0.485 mb, with NLO adding 0.011 mb, and N 2 LO another 0.007 mb. Cut-off dependence 
is negligible. The typical size of the expansion parameter in the pion-less EFT is about 
7 < /m 7r w 1/3. We therefore estimate the uncertainty from leaving out corrections at N 3 LO 
and higher as about 1/3 of the N 2 LO correction or 0.003 mb. 

IV. CONCLUSION 

The cross section for radiative capture of neutrons by deuterons nd — > 7 3 iJ at zero 
energies with was calculated pionless Effective Field Theory, the unique, model independent 
and systematic low-energy version of QCD for processes involving momenta below the pion 
mass. We applied pionless EFT to find numerical results for the Mi contributions. Incident 
thermal neutron energies have been considered for this capture process. At these energy 
our calculation is dominated by only S-wave state and magnetic transition M\ contribution. 
The M\ amplitude is calculated up to Next-to-Next to leading order N 2 LO. Three-Nucleon 
forces are needed up to N 2 LO order for cut-off independent results. The triton binding 
energy and nd scattering length in the triton channel have been used to fix them. Hence the 
cross-section is in total determined as a tot = [0.485(LO) + O.Oll(NLO) + 0.007(N 2 LO)] = 
[0.503 ± 0.003]m6. It converges order by order in low energy expansion. It is also cut-off 
independent at this order. We notice that our calculation has a systematic uncertainty from 
higher-order terms which is now smaller than the experimental error-bar. 
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